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Abstract 

Brown, Preston, and Singleton (BPS) produced an analytic calculation for energy exchange pro- 
cesses for a weakly to moderately coupled plasma: the electron-ion temperature equilibration rate 
and the charged particle stopping power. These precise calculations are accurate to leading and 
next-to-leading order in the plasma coupling parameter, and to all orders for two-body quantum 
scattering within the plasma. Classical molecular dynamics can provide another approach that can 
be rigorously implemented. It is therefore useful to compare the predictions from these two meth- 
ods, particularly since the former is theoretically based and the latter numerically. An agreement 
would provide both confidence in our theoretical machinery and in the reliability of the computer 
simulations. The comparisons can be made cleanly in the purely classical regime, thereby avoiding 
the arbitrariness associated with constructing effective potentials to mock up quantum effects. We 
present here the classical limit of the general result for the temperature equilibration rate presented 
in BPS. In particular, we examine the validity of the ?7ieiectron/"iion limit used in BPS to obtain 
a very simple analytic evaluation of the long-distance, collective effects in the background plasma. 
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Brown, Preston, and Singleton [ij] (BPS) have performed a controlled first principles 
analytic calculation of Coulomb energy exchange processes for weakly to moderately coupled 
plasmas, namely the electron-ion temperature equilibration rate and the charged particle 
stopping power. These calculations are accurate to leading and next-to-leading order in the 
plasma coupling parameter g ~ (e^/ADebye)/^, and to all orders in the quantum parameter 
fj ~ e^/hv-Y for the two-body quantum scattering within the plasma, where vt is a thermal 
averaged electron velocity. Here we shall examine the electron-ion temperature equilibration 



0,a. 



rate in some detail, particularly with regard to recent molecular dynamics calculations 

In general, we write the energy density exchange rate between the electrons and ions in 
a plasma as ^ 

^ = -a,(Te-TO, (1) 

where the rate parameter that appears here, 

/^y/' (2) 

2tx \2nTeJ ' ^ ^ 

serves to define the "Coulomb logarithm" In A for this specific process. The electron Debye 
wave number = l/X'^^ehyT^ and 

is sum of the squares of the ion plasma frequencies. Throughout this paper we shall measure 
temperature in energy units. 

For plasma conditions that often occur in ICF capsules, the coupling g is weak to mod- 
erate and the electron temperature Tg is much greater than the binding energy en of the 
hydrogen atom, Tg ^ Ch. Since v-^ ~ ^/Tjm'e, the quantum parameter can be alternatively 



be written as Tf ~ ^^ tYi/T^. Thus, for these ICF conditions, the quantum parameter is small, 
7] <^ 1, which corresponds formally to h being large, and so the scattering must be described 



^ We use the rate of change of electron-ion energy density d£/dt on the left hand-side of Eq. ([T]), rather 
than the corresponding temperatures dT/dt, because, in general, the conversion between £ and T entails 
higher order plasma corrections in the equation of state. However, as noted in Rcf. 1], a clean separation 
into separate electron and ion energy components within the plasma can be made only up to the order g 
to which we work here. The kinetic energy does not depend upon the coupling and its time derivative 
is of order In (7. On the other hand, the potential energy starts at order , and its time derivative is 
of higher order than we are calculating. At higher orders, collective, plasma potential energies come into 
play, and this clean separation into electron and ion energies cannot be made. 
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quantum mechanically. The electron mass me is much smaller than an ion mass rrii in the 
plasma. Taking advantage of rrie/mi <^ 1, BPS made use of a sum rule for the contribu- 
tions of collective long-distance effects in the me — >■ limit and were able to find a simple 
expression for the Coulomb logarithm under these conditions: 



HS^-7-1|, (4) 



where 7 = 0.57721 ■ ■ ■ is the Euler constant, and Ue is the electron plasma frequency. The 
result (jl]) corresponds to Eq's. (3.61) and (12.12) of BPS l|, with a small transcription error 

Cq. (12.43) to Eq. (12.44). 

have been used to extend such perturbative results 
into non-perturbative regimes. However, most MD methods employ ad-hoc potentials that 
mock up quantum-mechanical effects in a fashion that has little theoretical basis. Con- 
sequently, even in the perturbative regime, it is problematic to compare such MD results 



corrected when BPS passed from their '. 

r 

Molecular dynamic (MD) methods [2 



with the result of a systematic ca^ 
problem, Dimonte and Daligault 



culation such as that given in Eq. (jl]). To surmount this 
3| (DD) have recently performed purely classical ^ MD 
simulations with the Coulomb potential. The MD simulations of DD give results that can 
be directly compared with the classical limit of the BPS calculation. As with the extreme 
quantum limit (jll), for light electrons the classical result also takes a simple form. 

This expression follows from adding terms (12.25) and (12.44) of BPS, except for the trivial 
modification to ordinary Gaussian units from the rationalized electrostatic units employed 
by Ref. py]. Here, e and Cj = e are the values of the electron and ion charges. For 
the electron-proton plasma considered by DD, the classical BPS result can be written as 
InAgpg = ln{C/(7e}, where ge = e^Kg/Te is the plasma coupling, and C = 4e^*^^^+^/^) = 
0.7648 ■ ■ ■ . For small values of Qe, the regime in which the BPS formalism is valid, DD have 
found agreement with expression to within their statistical error of 5%. 

To avoid the formation of unstable configurations of two particles of opposite charge 
that experience an attractive force, DD work with electrons and ions of the same charge 



^ Dimonte and Daligault have obtained unambiguous and accurate results for a well-defined problem. They 
do not treat quantum-mechanical effects, and thus their work is not directly applicable for the conditions 
that appear, for example, in an ICF capsule. 



3 



moving in an implicit neutralizing background charge density of opposite sign. Since it is 
either the square or the absolute value of the electron and ion charges that enter into the 
perturbative regime, as made explicit by the expressions ([2D and this sign change causes 
no problem for small coupling ^. To increase the speed of convergence of the simulations, 
DD use an 'electron' mass rrig and a single species of 'ion' mass mj whose ratio is on the 
order of 1/100 rather than the physical value of 1/1836 for the electron and proton. This 
larger choice of electron-ion mass ratio allows them to obtain accurate numerical results 
with high statistics, and it in no way alters the relevance of their work in testing analytic 
results for which this mass ratio can be changed at will. However, to properly compare with 
BPS in this mass regime, one must now include finite electron mass corrections to Eq. (jS]). 
Some are trivial kinematic corrections. Others, which we shall denote by A, arise from long 
distance collective effects and are non-trivial. For the parameter regime considered in DD 
jsl, we shall find that these corrections are less than their 5% statistical error. 

Before turning to the calculation of the electron mass corrections for the classical case, we 
note that the extreme quantum limit (jlj) contains exactly the same non-trivial long distance 
correction A, in addition to other trivial kinematic corrections. Therefore, the calculation 
in the text is also relevant for the extreme quantum case. For plasmas in the IGF regime, 
these corrections may be of comparable size to the degeneracy corrections ^ calculated in 
Ref. 5|. To have a place in which all the small corrections to the quantum limit (jl]) are 



collected, we quote results of Ref. [5] in the Appendix. These entail not only the first 
Fermi-Dirac correction when the electron density starts to become large, but also the first 
classical correction when the quantum parameter rj starts to become large. Thus, we shall 
have a convenient reference to all of our results on the energy relaxation rate relevant to 
IGF recorded in the Appendix. 

To begin our development, we note that the work of BPS [l| expresses the coefficient 



^ However, in the non-perturbative regime where the plasma coupling g becomes of order one or larger, 
the relative sign of the electron and ion charges is important. For example, the first correction to the 
perturbative results quoted in the text has an overall coefficient involving (eei)^ that is negative in the 
physical case. This is experimentally confirmed in the differing ranges of positive and negative pi mesons 
stopping in nuclear emulsions Q| . 

Section III of Ref. [Sj contains a detailed pedagogical account of the method of dimensional continuation 
that we employ. Section IV of this paper reviews the ad-hoc nature of the "Convergent Kinetic Equations" 
that have been used in the literature. The work of Kihara and Aono Q which examines the energy 
relaxation of a slow particle in a plasma falls into this latter category. 
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appearing in Eq. ([T]) as the sum of three contributions, 

Ce:=(cr:,s + C,^:,H)+0, (6) 

with the three terms given by Eq's (12.25), (12.31), and (12.50) in BPS. The first term C^j g 
is the short- distance classical contribution. The second term C^^ ^ is the contribution from 
long-distance, cooperative, dielectric corrections. These two terms comprise the complete 
classical rate coefficient, which is the main topic of this note. The remaining term C^^ is the 
complete quantum correction (which vanishes in the formal limit h 0). The method used 
by BPS enables these terms to be unambiguously and precisely calculated to leading and 
first subleading order in the plasma coupling — expressed in terms of the plasma number 
density n, this is a unique evaluation to the formal orders n In n and n. Since the expression 
for C^j*^ is not needed for the discussion in the text, and since its ingredients are somewhat 
complex, we shall relegate it to the Appendix. The quantum result is, however, an essential 
contribution in practical ICF applications. 

The first term in Eq. ([6]), the short-distance classical scattering contribution, reads ^ 



(7) 



Here, the sum runs over all ions of charge Zi in the plasma, Qe is the electron coupling 
parameter, which in ordinary cgs units (not the rationalized units employed by BPS) is 

9e="4^. (8) 
e 

The square of the thermal velocity entering Eq. ([7]) is defined by 

T T 

= — + — , (9) 

rrii 



while the quantity 



(10) 



nXpi rrip rrii 



is the reciprocal of the reduced mass of an electron-ion pair. 



^ This is Eq. (12.25) of Ref. 1] in somewhat different notation, and with the arbitrary wave number K set 
to the electron Debye wave number Kg for convenience; see also Eq. (B8) of Ref. 
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The second term in Eq. ([6]), the long-distance, dielectric term that accounts for collective 
effects in the plasma, is given by 



pc ^ n ^ Pe{v)pi{v) I j i^totai(t^) , 

The individual spectral weights in this expression are defined by 



Pb{u) = kI ^^^^u exp{-mbU^/2Tb} , (12) 
the sum over ion components is denoted by 

Pi(w) = y'.P*(M) , (13) 
and so the total spectral weight is given by 

Ptotal(M) = Pe{u) + Pi{u) , (14) 

where Pe{u) is the electron contribution. The function -Ftotai('i^) is related to the classical, 
leading-order plasma dielectric permittivity by 

kh{k,kv) = e + FtoUv), (15) 

and it can be written in terms of the dispersion relation 

F,oUv)= r du ^^°^-'(^) (16) 

As explained in BPS [l| in the work about their Eq. (12.44), in the limit of small electron 
mass the long-distance contribution reduces, in leading order in mg, to the simple form 

.2 / _ \ 1/2 



Hl + ^l^^ (17) 



in which F^iy) differs form Ftot&iiy) by simply replacing ptotai(M) in the dispersion relation (fT6|) 
by Pi{u). The integrand in the leading term f|T7j) is analytic in the upper-half plane and 
hence, noting the asymptotic behavior -Fi(f) — * —Loljv^ as f — > oo, the integral can be 
evaluated by contour integration with the result that 

d L = -1 ^ f ^2 (IgN 

^^'^ 2 27r V 27r rj ' ^ ' 
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This is the result (12.44) of BPS after correcting for the trivial transcription error noted 
above. We express the general long distance contribution ffTTl) in the form 

CSh = ^ f^V'^f |-- + a|, (19) 
^^'^ 27r V 27r tJ M 2 ) ' ^ ^ 



which, upon comparing Eq. ( ITTl) with Eqs. ( ITTl) and ( JTSl) . defines 
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u^t J-oo 27r 



Ptotai(^^) fi:^ I f^l 



k2 



(20) 



The correction A vanishes as mg — ^ 0. All the various ratios that appear within the outer 
square brackets are dimensionless. The prefactor i^l/uj'^ has the dimensions an inverse veloc- 
ity squared, which combines with the integration measure dvv to produce a dimensionless 
quantity. Hence A is dimensionless as it must be, and therefore it is a function only of 
dimensionless quantities. One might expect that the only relevant dimensionless parameter 
is the ratio of the squares of the thermal velocity of the ions and electrons, Time/Terrii. 
However, the Debye wavenumbers are important and they involve the temperature and the 
density. Therefore, with the ion species at a common temperature Ti, the most general set of 
dimensionless quantities is the ion-electron temperature ratio Ti/T^, the ion-electron number 
density ratios ni/ue, the ion-electron mass ratios rrii/me, and the dimensionless ionic charges 
Zi = Cj/e. Hence, 

A = a(I^,1'M,1^].{Z,}). (21) 



T I I 

In the case of a single ion species, as considered by Dimonte and Daligault in Ref. |3|, 
the correction A depends only upon Tj/Tg, mi/mg, and Z^. We have examined analytic 
approximations for A, but they are long and cumbersome and do not provide insight into 
its structure. Hence in what follows, we shall present graphs of A obtained numerically for 
various parameters. 

Adding Eq's. ([7]) and ( flQll . and comparing with the definition (j2]) of the Coulomb loga- 
rithm shows that in the classical limit 



^ Note that the resuh of Ref. [gJ contains nothing in the nature of the finite electron mass, long-distance 
collective plasma contribution A. 
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FIG. 1: (Color) The BPS Coulomb logarithm In A^pg plotted as a function of me/rrii for four values 
of the electron-ion temperature ratio, Te/Tj = 0.5 ,1,5, 10, all with the coupling = 0.1. 

This reduces to Eq. ([5]) in the limit mg — 0. Dimonte and Daligault 3| use the conventional 
definition of the Coulomb logarithm for a single ion species rather than convention ([21), one 
that applies to plasmas containing a variety of ions. For a single species of ions, the two 
conventions are related by 

TiUle 



InA^ 



1 + 



-3/2 



In AC 



(23) 



where we have denoted the conventional definition'' of the Coulomb logarithm by InA*^. 
Pulling together previous definitions gives 



lnACp3 



In 



+ In 



rrij + rrie 



1 + 



Time 
Temi 



- 27 + 1 + 



TelTLi 



3/2 



1 



A 



(24) 



which in Fig. [T] is plotted as a function of the electron- ion mass ratio for several values of the 
temperature ratio. Upon expanding to leading order in mg/mi we can express the Coulomb 
logarithm in terms of a zero electron-mass contribution and a correction, 

(25) 



InAgps - 



lnA;^^° + AlnACp3 , 



^ As is apparent from Eq. (j22p . the logarithm depends upon ion species, and thus an overall factor of the 
form (1 + Time/TeTrii)^^^^ cannot be extracted in the general case. 
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FIG. 2: (Color) The correction A defined in Eq. (j20p plotted as a function of the mass ratio mg/mi 
for a hydrogen plasma. The four curves correspond to the temperature ratios T^/Ti = 0.5, 1, 5, 10. 



where 



lnA^P° = H^}-27-^ (26) 



is the zero electron- mass limit, and 



is the leading order electron mass correction. Dimonte and Daligault [3j use = 0.1 and 
consider the cases in which Tg/Tj varies from about 1 to 5 with me/rrii varying from about 
zero to 0.02, while Zi = 1 = Ue/rii are fixed. Figure [2] displays the values of A about 
this parameter range. Figure [3] presents the complete leading corrections for the Coulomb 
logarithm (|2^ as the mass ratio rrie/mi is varied — Eq. (!27j) divided by Eq. ( |26l) . The 
leading term (126|) is about 2.0, and so the relative correction is correspondingly smaller. 
Figure [3] shows that the relative size of the electron mass correction in the range examined 



by Dimonte and Daligault {s] is less than 2%, which is less than their statistical accuracy 
of 5%. With smaller statistical error, one could resolve the mass effects (127|) with an MD 
simulation. 
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FIG. 3: (Color) The relative electron-mass correction A In Agpg / In A^pg plotted as a function of 
the mass ratio me/rrii for four values of the temperature ratio T^/Ti = 0.5,1,5,10. In each case 
Qe = 0.1. The mass correction AlnAgpg is defined in Eq. (127]) . while the zero- mass logarithm 
In Agpg is given by Eq. ([26j) . 
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APPENDIX A: RESULTS COLLECTED 



We put together here all the formulae relevant to the cases of physical interest in which 
the scattering is predominantly quantum mechanical. These include the small classical 
corrections to this purely quantum limit and the leading effects of Fermi-Dirac statistics 
that come into play as the electron density is increased. This we do because, with the 
inclusion of the A correction, we now have in hand all the small corrections to the leading 
quantum-mechanical scattering limit. To exhibit these, we write 

In Abps = In A;^^3 + In A^^^ + In Al^^ . (Al) 

Here InAj^^g is the leading term in the quantum limit together with the A correction that 
we have exhibited in the text, InA^pg is the first classical correction that appears when the 
parameters depart from the extreme quantum limit, and InAgpg is the first correction when 
Fermi-Dirac statistics start to become important. The latter two terms have been computed 
in Ref. [sl to leading order in the small ratio rrie/mi] this suffices since the terms are already 
themselves small. 

In the text we examined the limit of purely classical scattering and thus omitted the 
quantum correction term C^^ in Eq. ^ from the complete relaxation rate. As a first step in 
presenting the collection mentioned, we quote this omitted correction which is Eq. (12.50) 
of BPS: 

3/2 



(A2) 



Here ip{z) = dlnr{z)/dz, and 



= m (A3' 



ei 



makes precise the definition of the quantum parameter alluded to at the beginning of the 
text with the square of the thermal velocity in this expression V"J defined previously in 
Eq. ([9]). The extreme quantum limit in which f/ej ^ of this formula is spelled out in detail 
in Ref. Here we shall not repeat the derivation but simply quote the BPS limit (12.53) 
with slightly different notation: 

^" = >^.^(r,.n,+T.m.)3/^l^J r + ^n^/J ■ ^^'^ 
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This quantum correction, added to the classical scattering contribution ([7]), produces 



C. 



el,S 



1/2 



3/2 r 



In 
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(A5) 



2 {Tme + T.mif V2vr. 
In the same way that the classical Coulomb logarithm (!22l) was constructed, the quantum 
scattering version now reads 



InA"^^ = V — I 1 + 



Time 
Terrii 



-3/2 



The explicit electron-ions mass ratio terms that appear here (including those contained in 
the definition of and Vj ) are easy to compute. For typical ICF conditions, they make 
very small corrections on the order or less than 1%. So as to make the significance of the A 
correction clear, a correction that does require some computation, we now omit these small 
explicit terms and write 



ill 1 Igpg 2 



In 



-7-1 



+ A 



(A7) 



which is precisely the result (jl]) of the text, but with the additional finite electron mass 
correction A. 

We show the A correction in Fig. H] over a wide range of the temperature ratio Tj/Te 
for the typical ICF case of an equimolar DT plasma. For a burning plasma, the Coulomb 
logarithm has the rough value In Aj^^g 4, and so the relative A correction is about a quarter 
of the number shown in Fig. HI 

For the remaining terms in Eq. ( lAll) . we shall just quote the result from Eq. (2.6) pre- 
sented in Ref. jsl, namely 



and 



In A™ 



InA^ 



rieK 



_ £h \ ^ 



2 , ,2 r 



- 1 



1 



23/2 



C(3) In 



-7 -2C'(3) 



In 



7-1 



1 , 1 
- In 2 + — - 

2 2^/2 



(A8) 



(A9) 



The ratio e-n/Te describes the relative size of the first quantum to classical correction, where 

13.6 eV (AlO) 



4 

e rrie 



2t9 



^ Reference fsj also contains the result in which the quantum-mechanical scattering is computed with exact 
Fermi-Dirac statistics, not just the first correction away from Maxwell-Boltzmann statistics which we 
quote here. 
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FIG. 4: (Color) The correction A defined in Eq. (j20p for an equimolar DT plasma (solid) and a 
hydrogen plasma (dashed) plotted as a function Tg/Ti for the physical values of the electron and 
ion masses. 



is the binding energy of the hydrogen atom. The numerical values of the zeta-function and 
its derivative are 

oo ^ 

C(3) = 5^ = 1.20205- ■■, (All) 



k=l 



and 



oo ^ 

('(3) = — lnA; = -0.19812--- . (A12) 
The electron thermal wave length 



k=l 



X, = h{^ (A13) 

\meTJ 

sets the scale at which quantum statistics comes into play, with Ze = neXl/2 the electron 
fugacity. 

For some temperature and number density regimes of interest, the two corrections flASp 



13 



and (1A9P become comparable in size 
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